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ABSTRACT 

The globular cluster co Centauri is one of the largest and most massive members of the galactic 
system. However, its classification as a globular cluster has been challenged making it a candidate 
for being the stripped core of an accreted dwarf galaxy; this together with the fact that it has one 
of the largest velocity dispersions for star clusters in our galaxy makes it an interesting candidate for 
harboring an intermediate mass black hole. We measure the surface brightness profile from integrated 
light on an HST/ ACS image of the center, and find a central power-law cusp of logarithmic slope 
-0.08. We also analyze Gemini GMOS-IFU kinematic data for a 5x5'' field centered on the nucleus 
of the cluster, as well as for a field 14'' away. We detect a clear rise in the velocity dispersion from 
18.6 km s-^ at 14" to 23 km s ^ in the center. A rise in the velocity dispersion could be due to 
a central black hole, a central concentration of stellar remnants, or a central orbital structure that 
is radially biased. We discuss each of these possibilities. An isotropic, spherical dynamical model 
implies a black hole mass of 4.01^ 0^ ^ lO^M©, and excludes the no black hole case at greater than 
99% significance. We have also run fiattened, orbit-based models and find similar results. While our 
preferred model is the existence of a central black hole, detailed numerical simulations are required to 
confidently rule out the other possibilities. 

Subject headings: globular clusters: individual (cc; Centauri), stellar dynamics, black hole physics 



1. INTRODUCTION 

The globular cluster co Centauri (NGC 5139) is re- 
garded to be the largest and most massive member of 
the Galactic cluster system with a tidal radius of 69 par- 
sees (lHarrisl [l996). an estimated mass of 5.1 x lO^M©, 
and a m easured central velo city dispersion of 22 ± 4 
km s~^ (|Mevlan et al.l [19951 ). The cluster presents a 
large scale global rotation, measured with radial veloc- 
ities, of 8 km s"-*^ a t a radius of llpc from the center 
( Merritt , et al.lll997 ) and c onfirmed with proper motions 
( van Leeuwen et aLl l2000l ). which make s it one of the 
most fiattened galactic globular clusters ([White fc Shawll 
[1987). A rotating fiattened model incl uding proper mo- 
tion and radial velocity datasets by Ivan de Ven et al.l 
(2OO6) calculate a lower total mass of 2.5 x IO^Mq and 
confirm the central line of sight velocity dispersion value 
of 20 km s-^ They measure a dynamical distance of 
4.8 ±0.3 kpc (which we adopt for this pape r), uj Cen has 
a peculiar highly bound retrograde orbit ( Dinescu et al.l 
I2OOI). It also has a stellar population that makes it stand 
out from the rest of the Galactic globular clusters due to 
its c omplexity. It shows a broad metall icity distribu- 
tion (jBedin et al.ll2004l : iNorris et al"]ll996D . as weh as a 
kinematical and spatial separation between the different 
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subpopulations ([Pancino et al.ll2003l : [Norris et al.lll997D . 

All the above results have led to the hypothesis that 
CO Cen is not a classical globular clust er, but instead 
is the nucleus of an a ccreted g:alaxv ([jFreemad Il993[ : 
iBekki fc FreemanI 120031 : iMeza et al.l 120051 ). The scenar- 
ios of it bein^ the product a merger of two globular 
clusters (Icke & Alcaino 1988) and of self-enrichment 
(Ikuta & Arimoto 2000) have also been proposed to ex- 
plain the stellar populations. 

The high measured velocity dispersion together with 
the possibility of being a stripped galaxy make uj Cen 
an interesting candidate for harboring a black hole in its 
cent er. The extrapo lation of M, — a relation for galax- 
ies ([Gebhardt et all l2000a[ : iFerrarese fc Merritt[ I2OOO : 
iTremaine et al.ll2002[ ) predicts a 1.3 x 10^ Mq black hole 
for this cluster. The sphere of infiuence of such black 
hole for a star cluster at the distance of lj Cen with a 
velocity dispersion of 20 km s~^ is ~ 5'', making it an 
excellent target for ground-based observations. 

Two globular clusters have been suggested for har- 
boring an intermediate mass black h ole in their nu- 
cleus. One is the g alactic cluster Ml 5 (Gebhar dt et al.l 
l2000b[ : iGerssen et a l. 2002, 2003) and t he other is Gl, 
a_giant^globular cluster around M31 (1 Gebhardt et al.l 
I2OO2L [2OO5I). M15 is assumed to be in a post-core 
collapse state, therefore its dynamical state has been 
debated between harboring a black hole or contain- 
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Fig. 1. — Surface brightness profiles for uj Cen. The circles 
show our measured photometric points from the ACS (filled) and 
H-alpha (open) images. The triangles show photometric points 
obtained from ground based images by Trager et al. The dashed 
line is Trager's Chebychev fit. The solid line is our smooth fit to the 
combination of the ACS points inside 40^^ and Trager's Chebychev 
fit outside. 

ing a lar^e number of comp act remnants in its cen- 
ter (|Baumgardt et aP l2QQ3bl lal) . Unfortunately, obser- 
vational constraint s between these two hypo thesis re- 
main inconclusive ( Ivan den Bosch erall l2QQ6l Gl on 
the other hand, has a core with characteristics closer 
to those of uo Cen, and obser vations support t h e black 
hole interpretation for Gl. iBaumg ardt et alJ (|2QQ3c[ ) 
propose an alternative interpretation for Gl in which 
they match the observations with a model of two col- 
liding globular clusters. The Gl black hole models are 
preferred since the M/L profile is expected to be flat 
in its core, so any rise in the velocity dispersion is un- 
likely to be due to remnants that concentrated there from 
mass segregation. The black hole interpretation for Gl 
is strongly supported by radio (Ulvestad et al. 2007) and 
x-ray (|Poolev fc Rappaportll2006l) detections centered on 
the nucleus, uo Cen and Gl have similar properties in 
both their photometric and kinematic profiles. In this 
paper we report photometric and kinematical measure- 
ments that suggest the presence of a central black hole 
in uo Cen. 

2. SURFACE BRIGHTNESS PROFILE 

The surface density profile at large radii for uo Cen has 
been measured from a combination of star count s and 
apert ure photometry froni ground based images (.Mevlanl 
119871 : iTrager et allllQQa Ivan Leeuwen et al.N2000f ). We 
measure the central part of the profile taking advan- 
tage of iJ-i/ 66 Space Telescope (HST) spatial resolution. 
We measure integrated light from an ACS F435W im- 
ag e (340 sec) applying the te chnique described in detail 
in iNovola fc Gebhardtl (|2006l ). which uses a robust sta- 
tistical estimator, the bi- weight, to calculate number of 
counts per pixel on a given annulus around the center 
of the cluster. As a test, we also measure the profile 
from a narro w band H-alph a image from the Rutgers 
Fabry-Perot (jXie et al.ll2006l ) with lower spatial resolu- 



tion. Since both images ha ve a limited r adial coverage, 
we use the Chebychev fit of iTrager et al.l (pJ95) for the 
surface brightness profile to cover the full, radial extent 
of the cluster. All profiles are normalized to the Trager 
profile. 

Having accurate coordinates for the center of the clus- 
ter is crucial when measuring density profiles. Using the 
wrong center typically produces a shallower inner profile. 
We use a technique where we take an initial guess cen- 
ter, divide the cluster in eight concentric sectors around 
this center, and calculate the standard deviation of the 
sum of stars for the eight sectors. The radius of the sec- 
tors is chosen to be as large as the image will allow, in 
this case it is ~2'. We repeat the procedure for a grid 
of center coordinates and use the one that has the min- 
imum standard deviation. Details a bout the technique 
can be found in iNovola fc Gebhardtl (j2006). The coor- 
dinates for our center are RA 13 : 26 : 46.043 and DEC 
-47 : 28 : 44.8 on the ACS dataset J6LP05WEQ using 
its WCS zeropoint. 

The measured profiles from the B-band and H-alpha 
images are consistent as can be seen on Figure [H The 
H-alpha profile follows the turnover around the core ra- 
dius very well up to 100^' and it also shows the rise toward 
the center, but it is noisier than the ACS profile. The 
solid line is a smooth fit made to the combination of the 
photometric points from ACS inside 40'' and Trager et al. 
Ch ebychev fit outside 4 0''. For comparison, we include 
the ITrager et al.l (|l995l ) photometric points in the plot. 
The surface brightness profile shows a continuous rise to- 
ward the center with a logarithmic slope of —0.08 ±0.03, 
which is in contrast to the common notion that uo Cen has 
a flat core. Van Leeuwen et al (2000) perform star counts 
for giant stars and notice that they are more concentrated 
than previously thoug ht . Our r esult is consistent with 
their finding. Baumgar dt et al.l (|2005l ) perform N-body 
models of star clusters with initial King profiles and con- 
taining a central black hole. They predict the formation 
of a shallow cusp of —0.1 t o —0.3 logarithmic slope af- 
ter 1.5 — 4 relaxation times. Trenti et al. (2007) perform 
similar N-body models and conclude that clusters con- 
taining central IMBHs evolve to configurations that have 
extended cores, with rc/vh ~ 0.3. The rc/vh ratio for uo 
Centauri is 0.33, which is in very good agreement with 
the predictions from this models. Our observed surface 
brightness profile is intriguing considering that it follows 
the two specific predictions from these N-body models, 
but of course, only kinematical measurements can deter- 
mine the mass profile, and test the existence of a central 
black hole. 

3. KINEMATIC MEASUREMENTS 

Obtaining kinematic information of the central regions 
of Galactic globular clusters is a challenging task since 
the brightest stars dominate with typical ground-based 
conditions and the extreme crowding produces confusion. 
Measuring individual radial velocities requires a spatial 
resolution that can only be achieved with adaptive op- 
tics or from space. On the other hand, measuring velocity 
dispersion from an integrated spectrum is subject to shot 
noise due to giant stars whose contribution dominates 
the light. D ubath et al. (1997 ) calculate the relative con- 
tribution to integrated light by different stellar groups. 
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Fig. 2. The central field observed with the GMOS-IFU. The 
green box represents the GMOS field of view. Top left: ACS image 
of the observed region. The red circle marks the center of the clus- 
ter. Top right: Convolved ACS image to reproduce the reported 
seeing during observations. Bottom left: GMOS acquisition image. 
Bottom right: Reconstructed GMOS-IFU image. 
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Fig. 3. — Same as fig 2 for the field 14'^ away. 



They find that the contribution from the few brightest 
stars is equal in weight to that of the much larger num- 
bers of fainter stars. Therefore, the only way to obtain 
accurate radial velocity dispersion measurements from an 
integrated spectrum is if the participation of the bright- 
est stars can somehow be avoided or at least minimized. 
One way to do this is by observing crowded regions with 
an integral field unit (IFU) which produces individual 
spectra of subsections in the region (typically ~0.2'' in 
size). One can exclude the spectra affected by the bright- 
est stars when measuring the integrated background light 
and thus decrease the shot-noise contribution to the un- 
certainty. The Gemini telescopes operate primarily us- 
ing a queue scheduling, which makes them an excellent 
tool to measure integral field spectra of globular clusters 
since observing constraints (such as excellent seeing) can 
be specified in advance, and data is only taken when the 
required observing constraints are met. 

As part of the Science Verification program for the 
Gemini GMOS-South IFU, we obtained nod- and- shuffle 
observations on February 29 2004 (program ID: GS- 
2003B-SV-208). We use the IFU in 2-slit nod- and- shuffle 
mode, which gives a field of view of 5'' x 5'', comprised of 
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Fig. 4. — Velocity dispersion profile for u; Cen with various central 
black hole models. Filled squares are the dispersions and uncer- 
tainties from the GMOS-IFU and open circles are from individual 
radial velocity measurements. A set of isotropic spherical models 
of varying black hole masses is shown for comparison. The thick 
line is the no black hole model and the thin lines represent models 
with black holes as labeled. The dashed (red) line is the velocity 
dispersion profile for the best-fit orbit-based model (in Section 4). 

700 individual lenslets plus fiber elements, each of which 
covers approximately 0''.2 on the sky. We use the R600 
grating, yielding a resolving power R=5560, which we 
measure from the lamp spectral lines, along with the 
Calcium Triplet filter to give a wavelength coverage of 
7900-9300A. Three fields are observed, each for a to- 
tal integration time of 900 sec on source and 900 sec 
on sky. The observations are made using the nod-and- 
shuffle technique with 30 sec sub-integrations observed 
in a B-A-A-B pattern, where A is the on-source posi- 
tion and B is the sky position, located 498^' away. The 
nod- and- shuffle technique improves the sky subtraction, 
especially in the presence of CCD fringing, by sampling 
the object and sky on exactly the same CCD pixels, with 
exactly the same light path, on timescales comparable to 
those of the sky emission line variability. The first of the 
three fields is located at the cluster center, and the sec- 
ond field is centered W away. The third field appears 
to have been pointing somewhere else but, despite much 
efforts, we cannot determine exactly where the IFU ob- 
servations are pointed (they do not match anything in 
the acquisition image for this field). The reconstructed 
IFU image for this third field contains fewer stars and 
the PSF is obviously broader than for the other two. It 
is clear that the exposure was taken during much worse 
seeing conditions than the other two fields, so shot-noise 
effects are likely to be important; for this reasons we 
exclude the third field from further analysis. Using the 
standard tasks from the IRAF- GEMINI package we sky 
subtract, fiat-field, and extract the spectra for each fiber 
and apply a wavelength calibration. 

The standard flat subtraction does not remove all of 
the fringing pattern in the image. As a result, a con- 
stant number of counts have to be subtracted from the 
data frames before flattening in order to reduce fring- 
ing problems. Relative to the bias frame, the amount 
of additional counts from scattered light is about 8%. 
Even after this procedure, there is some residual fring- 
ing that can only be alleviated by combining individual 
flbers over the full fleld into one spectrum. To combine 
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individual fibers we first divide by the continuum. The 
intention is to de-weight the bright stars with respect to 
the fainter ones, which helps to lessen the problems due 
to shot noise. For the continuum fit, we run a boxcar 
of dimension 111 x 1 over the reduced image, and then 
divide the central pixel by the median of the pixels in 
the box. This procedure brings all spectra to the same 
continuum level. We then combine every individual fiber 
with the six adjacent ones, since this represents one see- 
ing disk for the observations 

Figure [2] shows the reconstructed image from the IFU 
fibers for the central frame and the acquisition image 
as well as the same region on the ACS image. We also 
show a convolved image (with the reported seeing for the 
observations) of the ACS frame. The same match is per- 
formed for the field 14'' away (Fig[3|). Both ACS fields 
contain ~ 100 resolved stars. We construct a luminosity 
function for the detected stars for each field and com- 
pare it to the luminosity function of the entire cluster 
core. The luminosity function is consistent between the 
two fields. The brightest stars detected in both fields are 
two magnitudes fainter than the brightest stars detected 
in the core of the cluster. This excludes the possibility of 
the integrated spectra being artificially broadened by the 
presence of more blue straggler stars in the central field 
compared to the field 14'' away. Using the photometric 
measurements of individual stars together with the re- 
ported seeing, we calculate how many stars contribute 
to each fiber. Excluding the fibers which are dominated 
by a single star we estimate that the integrated spec- 
trum of the background unresolved light represent about 
60 stars in both fields. 

We focus on the Ca triplet region (8450A-8700A) for 
our analysis. We measure the relative velocities between 
each fiber for the two fields and obtain velocity distribu- 
tions from the individual fiber velocities. We fit a Gaus- 
sian to the velocity distributions and observe that the 
one for the central field is clearly broader than for the 
one 14" away. The largest relative velocity between two 
fibers is 80 km s~^ for the central frame, and 60 km s~^ 
for the other one. Using the dispersion of the individual 
fiber velocities as a measure of the cluster velocity dis- 
persion will be biased. Since multiple stars, in general, 
provide light to an single fiber, the measured velocity in 
that fiber will be pulled toward the cluster mean as op- 
posed to representing one star. Thus, the dispersions of 
the fiber velocities will be biased significantly low. This 
is what we find although the central frame does have a 
obviously larger spread in fiber velocities. 

To properly estimate the velocity dispersion we have 
to rely on the integrated light, and require template stars 
in this case. Unfortunately, we do not have isolated stars 
that are free from the fringe problems mentioned ear- 
lier, so we cannot accurately use template stars observed 
with the sam e instrument. We rel y on the template stars 
observed bv Wal cher et all (|2005r ). from VLT-UVES ob- 
servations at around R=35000. We convolve the spectra 
to our measured resolving power. To extract the veloc- 
ity dispersion from the integrated light we utilize the 
non-parametric, pixel- based techniqu e as described in 
([Gebhardt et al. 20003 iPinknev et al.| !2003). We choose 
an initial velocity profile in bins, convolve it with the 
template (or set of templates), and calculate residuals to 
the integrated spectrum. The parameters for the line of 
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Fig. 5. — vs. black hole mass. The minimum is found for 
a black hole mass of 4.0 x 10^ Mq, with 68% confidence limit at 
3 and 4.75 x 10^ marked by the dashed line. For our model 
assumptions, the no black hole model is excluded at greater than 
the 99% confidence 

sight velocity distribution (either velocity bin values or, 
if desired, a parametric Gauss-Hermite expansion) are 
varied to provide the best match. Monte Carlo simula- 
tions determine the uncertainties, and use the measured 
noise in the spectrum. 

The dispersion fitting routine allows for a mismatch in 
the equivalent width between the object and template. In 
this case, there is a 30% difference in the equivalent width 
of the calcium triplet lines. We do not know whether 
this is caused by the scattered light (unlikely given the 
amplitude), the specific templates we used, omega Cen's 
particular composition, or a combination of all three. We 
have run a variety of tests to determine whether stars of 
difference equivalent widths would cause a bias in dis- 
persions, and find no such bias. We have also measured 
the dispersions using template stars from the same IFU 
data, since there is at least one star that is fairly iso- 
lated. The uncertainties are larger due to the scattered 
light problems, but the value of the dispersion remains 
the same. Thus, we conclude that template issues are not 
a significant source of bias in the dispersion estimate. 

We combine the spectra from individual fibers using a 
biweight estimator. Different sets of fibers for each frame 
are combined in order to test for consistency in our re- 
sults. First, we combine every fiber on the frame, then 
we exclude the 25%, 50% and 75% brightest fibers. We 
measure the velocity dispersion for these four spectra for 
each frame. The measured velocity for the central frame 
is always higher than the one for the frame 14" away 
for every equivalent pair of combined spectra. We mea- 
sure velocity dispersions from 21.8 to 25.2 km s~^ for 
the central field, and 18.2 to 19.1 km s"^ for the field 14" 
away. We adopt 23.0 ± 2.0km s~^ for the central field and 
18.6 ± 1.6km s~^ for the other. The latter measurement 
coincides with the central velocity dispersion value mea- 
sured for (jj Cen by various authors. Van den Ven et al. 
(2006) measure a line of sight velocity dispe rsion profile 
by combining vari ous datasets ( Suntzeff fc Kraft 1 19961 : 
iMavor et all 1 19971 : iReiins et al] l2006l : IXie et all l2006[ ). 
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Fig. 6. — Inferred density profiles. The solid line is the depro- 
jected density profile for the luminous component. The dashed line 
represents the required dark component to reproduce the observed 
kinematics. 

They use 2163 individual radial velocity measurements 
divided into polar apertures to obtain the final velocity 
dispersion profile. We use dispersion estimates as pre- 
sented in their table 4. The average radius for the spec- 
tra that contribute to the central values is 2.b" and lA" 
for the second. Figure [4] presents the velocity dispersion 
data. 

4. MODELS 

As discussed in Section 2, the central shape of the sur- 
face brightness profile of uo Cen resembles that found by 
[Baumgardt et al. (2005) in star clusters harboring black 
holes. The presence of an intermediate black hole at the 
center of this cluster is one of the possibilities for explain- 
ing the observed rise in velocity dispersion. We have run 
two types of modeling: 1) spherical and isotropic, and 
2) flattened, orbit-based models that allow for general 
anisotropy. Although the orbit-based models are more 
general, they do not consider dynamical stability. For 
a system with a half light relaxation time shorter than 
its age, like uo Cen, the dynamical evolution cannot be 
ignored. Therefore, while the two models give similar 
results we quote results from the isotropic analysis and 
use the orbit-based models to explore possible velocity 
anisotropics. 

For the isotropic analysis, we create a series of 
models using the no n-parametric method described in 
iGebhardt fc Fische3 |1995). As a first step we apply a 
reddeni ng correction t o the observed surface brightness 
profile. iHarrisI (|l996f ) reports a 0.1 reddening for this 
cluster, which, although being relatively low, it is impor- 
tant for the proper M/L determination of the models. 
The reddening correction will only affect the M/ L value 
of the models, but not the shape of the profiles. We 
deproject the surface brightness profile using Abel inte- 
grals assuming spheroidal symmetry. The integral in- 
volves a derivative of the profile, therefore, any amount 
of noise present is amplified. We apply a spline smoother 
to the surface brightness profile before deprojecting and 
thus obtain a luminosity density profile as discussed in 
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Fig. 7. — Contours of as a function of black hole mass and 
mass-to-light ratio. Each point represents a particular model. The 
contours represent the 68, 90, 95, and 99% confidence for one 
degree-of-freedom, implying Ax^ = 1.0,2.7,4.0, and 6.6. The cir- 
cled point is the model that has the minimum value. 



I Gebhardt enH (119961). By assuming an M/L ratio, we 
calculate a mass density profile, from which the potential 
and the velocity dispersion can be derived. We repeat the 
calculation adding a variety of central point masses rang- 
ing from to 7.5 x lO^M© while keeping the global M/L 
value fixed. Van de Ven et al (2006) measure a fairly con- 
stant stellar M/L profile for uo Cen of 2.5 ±0.1. We find a 
constant mass luminosity ratio of 2.7, since this provides 
the best match to the observed velocity dispersion profile 
outside the core. 

Figure [4] shows the comparison between the different 
models and the measured dispersion profile. The most 
relevant part of the comparison is the rise inside the core 
radius, in particular the rise between the two innermost 
measurements. As it can be seen, an isotropic model with 
no black hole present predicts a slight decline in velocity 
dispersion toward the center, instead we observe a clear 
rise. The predicted central velocity for the no black hole 
model is 14.6 km s~^ which is well below any line of sight 
velocity dispersion measured inside l^ The calculated 
values for each model are plotted in Figure [5l as well 
as a line showing a Ax^ = 1. The x 
a best-fitted black hole mass of 4+^ '^^ x 
the original velocity dispersion profile without our two 
innermost measurement already points to an intriguing 
discrepancy, but the central measurements confirm an 
important rise in M/ L from the core radius to the center 
of this cluster. The central M/L value is 6.7, which is 
a considerable rise from the value of 2.7 just inside the 
core radius. Our best fit model implies a central density 
of 5.6 X Mq/pc^ the largest measured in a globular 
cluster. 

We also construct axisymmetric orbit-based dynam- 
ical mo dels. The models a re based on the formula- 
tion bv |Sch warzschild' (1979*) and are constructed as in 
IGebhardt e t al. (2000c, 2003). These models provide the 
most freedom possible of the distribution function for an 
axisymmetric system. We use the same deprojection as 
describe above, except we also include the observed flat- 
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tening. Assuming an M/L ratio and a BH mass, the 
mass distribution of the cluster is obtained and from it, 
the potential can be computed. Using this potential, we 
generate about 10^ representative orbits. The best match 
to the observed photometric and kinematical data pro- 
vide the orbital weights for a given potential. The pro- 
cess is repeated for various M/L values and BH masses 
until the minimum model is found. The kinemati- 
cal observations matched here are the individual radial 
velocities from van den Ven et al. (2006), and the two in- 
tegrated measurements described above. We plan to per- 
form more detailed orbit-based models including proper 
motion measurements in the future. 

Figure [3 plots a function of black hole mass and 

stellar M/L. The best fit model requires a black hole 
mass of 3 ± 1 X 10^ M©, with 1-sigma of the isotropic 
result. As expected, since orbit-based models are more 
general, the difference in compared to the no black 
hole case is smaller than in the isotropic case. The Ax^ 
between the best-fitted black hole model and no black 
hole model is 4 (marginalized over M/L), implying a 95% 
significance. In comparison, the difference for the 
isotropic models is 25. As discussed below, the reason 
for the difference is that the orbit-based models allow 
for significant radial anisotropy in the core for the no 
black hole model, while the best-fitted mass produces 
a nearly isotro pic distribution consist ent with previous 
measurements (|van de Ven et al.ll2006[ ). 

We also plot the velocity dispersion profile as mea- 
sured in the best-fit orbit based model in Figured! given 
as the dashed (red) line. The dispersion profile for the 
orbit-based model with no black hole is very similar, with 
no obvious correlated differences. The reason there are 
no obvious differences is because the overall change in 
is small and the orbit-based models tend to redistribute 
orbital weights to spread x^ over the full radial range. 
As discussed in Gebhardt et al. (2003), it is difficult to 
see differences in radial profiles of projected kinematics 
between models with black holes and without, even for 
galaxies where the overall x^ difference is large. This is 
understood since the orbit models have the freedom to 
change the orbital properties in such a way to even out 
the x^ differences over the full spatial extent of the ob- 
servations. For omega Gen, the black hole model only 
provides an increase in x^ of 4; thus with 23 data points, 
the average difference in terms of the measurement uncer- 
tainty is 0.42. For these reasons, the orbit-based models 
provide a modest significance for a black hole, and the 
argument is significantly strengthened when considering 
the need for the strong amount of radial anisotropy when 
a black hole is not included. 

The observations we use in the dynamical modeling 
rely on only the first and second moments of the velocity 
distribution. Since we have individual velocities (except 
for the central two radial points), we can utilize the full 
velocity profile. Furthermore, proper motion data ex- 
ists for data at larger radii (van den Ven et al. 2006). 
Including both effects, full velocity profiles and proper 
motions, will be the subject of a more detailed paper on 
omega Gen. 

4.1. Alternative to a Black Hole: Dark Remnants 

A possible alternative to explain the observed rise in 
M/L toward the center is a concentration of dark stellar 



Fig. 8. — Radial over tangential anisotropy vs. radius from orbit- 
based models. The solid line is for the best fit model containing 
a black hole. The dashed line is for a model with no black hole 
present. The vertical line marks the location of the core radius. 

remnants, e.g. neutron stars, stellar mass black holes, 
or massive white dwarfs. Using the observed velocity 
dispersion profile, we calculate the total enclosed mass 
and from this, the mass density profile. We then com- 
pare this with the enclosed mass implied by the lumi- 
nosity density profile, assuming the same M/L as for the 
models in the previous section. From these two profiles, 
we can estimate the density profile of the implied ex- 
tended dark component, if we assume this was the cause 
of the velocity dispersion rise towards the center. Figure 
[6] shows the estimated density profile for the dark and lu- 
minous components. It is clear that, if the velocity rise is 
due to an extended distribution of dark stellar remnants, 
the density profile of this dark component needs to be 
very concentrated and steep, with a logarithmic slope of 
~ —2.0 (the slope of the dashed line in the figure), resem- 
bling a cluster undergoing core collapse. The relaxation 
time for uo Gen implies a much slower dynamical evolu- 
tion than the one necessary to reach such a configuration. 
Gore-collapse models have shown that when a cluster has 
reached such a high degree of mass segregation, the ob- 
servable core to ha lf light radius ^ets very small, wit h 
values below 0.05 (iBreeden et al.l Il994l : iMakinol fl996[ ) . 
while this ratio is 0.3 for uj Gen. Also, the concentration 
value for uo Gen is c = 1.6 (Harris 1996), which is too low 
a value for the cluster to have undergone core-collapse. 
No evolutionary model predicts such a concentrated dis- 
tribution of dark remnants inside a cluster with a shallow 
extended core for the visible stars. Furthermore, the re- 
quired number of dark remnants is around 1% of the 
cluster mass. While this number is expected from stellar 
evolution, it is not expected to have all of the remnants 
to be concentrated inside of the core. 

4.2. Alternative to a Black Hole: Orbital Anisotropy 

Another possibility is that the observed rise in velocity 
dispersion is due to velocity anisotropy in the cluster. 
The orbit-based models explore this possibility. The two 
main results from the orbit analysis are that 1) the lowest 
X^ model is the one with a central black hole and 2) the 
model with no black hole requires a substantial amount 
of radial anisotropy. Figure [8] shows that without the 
presence of a central black hole, a large degree of radial 
anisotropy - (7^/(7^=1.5 - is required inside 0.3 Tc. At r > 
28'', the models with and without a black hole are close 
to isotropic, in agreement with the results of van den Ven 
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et al. (2006). For a system as dense as uo Cen, such a 
degree of anisotropy as measured in the no black hole 
case is expected to be quickly erased through relaxation 
processes. However, even with such a strong amount of 
radial anisotropy, the no black hole case is a poorer fit 
than the best-fitted black hole. One of the disadvantages 
of the orbit-based models is that we cannot include the 
dynamical stability arguments in the analysis, and for 
a system with a short relaxation time such as a globular 
cluster, this may be important. We, therefore, adopt 
results from the isotropic models, in particular since the 
analysis of van den Ven et al. (2006) find an isotropic 
distribution. 

5. DISCUSSION 

We measure the surface brightness profile for the glob- 
ular cluster uo Centauri (NGC 5139) from an ACS image 
in the central 40^'. The profile shows a continuous rise to- 
ward the center with a logarithmic slope of —0.08 ±0.03, 
in contrast with previous measurements which found a 
flat core. The shape of the profile is similar to that ob- 
tained from numerical models of star clusters containing 
black holes in their centers. We measure a line of sight 
velocity dispersion for two xb" regions, one at the cen- 
ter of the cluster and the other lA" away. We detect a rise 
in velocity dispersion from 18.6 km s~^ for the outer field 
to 23 km s~^ for the central one. We combine these two 
measurements with previously measured velocity disper- 
sion at larger radii. 

When we compare the observed velocity dispersion 
profile with a series of isotropic models containing black 
holes of various masses, we find that a black hole of 
4.0t?;o^ X lO^M© is necessary to match the observa- 
tions. We explore alternative explanations for the ob- 
served rise in our central velocity dispersion measure- 
ments. First we consider the possibility that the observed 
M/L rise is due to the presence of an extended compo- 
nent composed of dark remnants such as neutron stars 
or faint white dwarfs. The density profile of the dark 
component is required to be extremely concentrated to- 
ward the center, with a configuration practically decou- 
pled from the luminous component, uo Cen has a weak 
cusp in the central luminosity density profile, implying 
that the gravitational potential is very shallow inside the 
core and ther efore mass segregation is only a weak effect. 
[F^rraro et al.l (|2006l ) confirm the lack of segregation by 
measuring the radial distribution of blue straggler stars, 
which are heavy stars and should sink to the center of the 
cluster if there is mass segregation. They find a flat ra- 
dial distribution of blue stragglers with respect to lighter 
stellar populations. Also, the formation channel of the 
blue stragglers is not collisional, as it would be expected 
if there was a considerable amount of mass segregation. 
With this evidence in hand, there is no reason to expect 
a large variation of M/L inside the core due to stellar 
content, so a detected rise in M/L is likely to come from 
the presence of a concentrated massive object. 

There is also a stability argument against a dense com- 
pact cluster of dark remnants. The central density as 
measured from the 23 km s~^ dispersion estimate at 1.8'' 
is 5.6 X l{f Mq/pc^ . This is the largest measured for a 
globular cluster and it would be difficult to maintain us- 
ing stellar remnants. Obviously, if the density is due 
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Fig. 9. — M» — dvei relatio n for elliptical galaxie s and bulges. 
The solid line is the relation in'Tremaine et al " ("2002). uo Cen lies 
on the low mass extrapolation and suggests a similarity between it 
and the galaxies. Different types of systems such as star clusters 
and low luminosity AGN appear to populate the low mass end of 
the diagram. 

to solar mass remnants, over 10^ remnants would be re- 
quired inside of 0.05 parsecs. Using the arguments of 
Maoz (1998) and Miller (2006), this mass and density 
makes uo Cen one of the better examples where stellar 
remnants can be ruled out due to evaporation. Maoz es- 
timates that for these numbers, any cluster of remnants 
will have evaporated within the age of the cluster. 

An observed velocity dispersion rise toward the cen- 
ter of a cluster can also occur if a degree of anisotropy is 
present. If more radial orbits are present, those stars pass 
near the center at higher velocities than they would in 
an isotropic case. This possibility is evaluated with our 
orbit-based models. Our results agree with the model by 
van den Ven et al. (2006) in showing no anisotropy in the 
central 10 arcmin. Their models show a degree of tangen- 
tial anisotropy at large radius, but no radial anisotropy. 
Models without a central black hole but having a large 
degree of anisotropy inside 28'' are not as as good as mod- 
els including a black hole. Furthermore, the high degree 
of radial anisotropy is highly unstable in dense systems 
like uo Cen and therefore it is an unlikely explanation for 
the observed kinematics. 

Figure [9] shows the known M, — cTy r elation for black 
holes i n elliptical galaxies and b ulges (|Gebhardt et al.l 
"2000a'; iFerrarese fc Merritt' '200 QD . The galaxies used to 
determine the relation (Tremai ne et al. ^ 2002) are plot- 
ted along with objects containing smaller black holes in 
low luminosity quasars (Ba rth et al. ,2005) . two nearby 
low luminosity AGN (NGC 4595 and Pox 52), and three 
globular clusters (Gl, Ml 5 and uo Cen). We also plot the 
upper limit for th e black hole mass in the nucleus of M33 
(|Gebhardt et al.ll2001l ). which does not lie on the correla- 
tion. The black hole in uo Cen lies above the relation, but 
it is consistent with the scatter observed a larger masses. 
The measured black hole mass is 1.6% of the total mass of 
the cluster, which is much larger than the canonical valu e 
of ~ 0.3% for larger spheroids (Magorrian et al.lll998f ). 
If uo Cen is indeed the nucleus of an accreted galaxy it is 
expected that it's original mass was considerably larger 
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than what we measure now. iBekki fc Fre^^^ ([2QQ3t ) 
reproduce the current mass and orbital characteristics of 
uo Cen with a model of an accreted lO^M© dwarf galaxy. 
A mass of 4 x lO^M© for the original spheroid would put 
the black hole near the 0.3% value. 

The two pieces of observational evidence that uo Cen 
could harbor a central black hole come from the photom- 
etry and the kinematics. From the HST image of uj Cen, 
we measure a central logarithmic surface brightness slope 
of —0.08 ±0.03. This value is very similar to that claimed 
by the N-body simulations of Baumgardt et al. (2005) 
that are most likely explained by a central black hole. 
Standard core-collapse does not lead to such a large core 
with a shallow central slope. The black hole tends to 
prevent core collapse while leaving an imprint of a shal- 
low cusp. It will be important to run models tailored to 
uj Cen to see if one can cause and maintain a shallow cusp 
without invoking a central black hole. However, the main 
observational evidence for the central mass comes from 
the increase in the central velocity dispersion, where we 
detect a rise from 18.6 to 23 km s~^ from radii of 14 to 
2.b" . In fact, even excluding the Gemini data presented 
here, the previous ground-based data suggest a central 
mass concentration as well. The core of uj Cen is around 
155 " (about 2.5 so the dispersion rise is seen well 
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